Phase space deformation of a trapped dipolar Fermi gas 
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We consider a system of quantum degenerate spin polarized fermions in a harmonic trap at zero 
temperature, interacting via dipole-dipole forces. We introduce a variational Wigner function to 
describe the deformation and compression of the Fermi gas in phase space and use it to examine 
the stability of the system. We emphasize the important roles played by the Fock exchange term of 
the dipolar interaction which results in a non-spherical Fermi surface. 
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Two-body collisions in usual ultracold atomic systems 
can be described by short-range interactions. The suc- 
cessful realization of chromium Bose-Einstein conden- 
sate (BEC) [H and recent progress in creating heteronu- 
clear polar molecules Q have stimulated great interest 
in quantum degenerate dipolar gases. The anisotropic 
and long-range nature of the dipolar interaction makes 
the dipolar systems different from non-dipolar ones in 
many qualitative ways [3|. Although most of the theo- 
retical studies of dipolar gases have been focused on dipo- 
lar BECs, where the stability and excitations of the sys- 
tem are investigated (see Ref. Q and references therein, 
and also Ref. [J] ) and new quantum phases are predicted 
[El @] ! some interesting works about dipolar Fermi gas do 
exist. These studies concern the ground state properties 
0,0, i], dipolar- induced superfluidity [ljj , and strongly 
correlated states in rotating dipolar Fermi gases 111. 
None of these studies, however, takes the Fock exchange 
term of dipolar interaction into proper account 12 1. 

In this Letter, we study a system of dipolar spin po- 
larized Fermi gas. We will show that the Fock exchange 
term that is neglected in previous studies plays a crucial 
role. In particular, it leads to the deformation of Fermi 
surface which controls the properties of fermionic sys- 
tems, and it affects the stability property of the system. 
As Fermi surface can be readily imaged using time-of- 
flight technique [li| , this property thus offers a straight- 
forward way of detecting dipolar effects in Fermi gases. 

In our work, we consider a trapped dipolar gas of single 
component fermions of mass m and magnetic or electric 
dipole moment d at zero temperature. The dipoles are 
assumed to be polarized along the z-axis. The system is 
described by the Hamiltonian 



H 



N 

E 

i=l 



"2^ + ^ 



2 V 

i¥>3 



vUn-rj), (i) 



where Fdrf(r) = (d 2 /r 3 )(l-3z 2 /r 2 ) is the two-body dipo- 
lar interaction and U (r) the trap potential. To character- 



ize the system, we use a semiclassical approach in which 
the one-body density matrix is given by 
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where f(r, k) is the Wigner distribution function. The 
density distributions in real and momentum space are 
then given respectively by 

n(r) = p(r,r) = (2tt)~ 3 J d 3 kf{r,k) , 

n(k) = (2ir)- 3 J d 3 rf(r,k) . 

Our goal is to examine n(r) and n(k), as well as the 
stability of the system by minimizing the energy func- 
tional using a variational method. Within the Thomas- 
Fermi-Dirac approximation 0], the total energy of the 
system is given by E — Eki n + E tr + Ed + E ex , where 
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The dipolar interaction induces two contributions: the 
Hartree direct energy E d and the Fock exchange energy 
E ex . The latter arises due to the requirement of the an- 
tisymmetrization of many-body fermion wave functions 
and is therefore absent for the dipolar BECs. 

Homogeneous case — Let us first consider a homoge- 
neous system of volume V with number density n / , which 
will provide some insights into the trapped system to be 
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FIG. 1: Deformation function I(x). 

discussed later. In this case, we obviously have E tr = 0. 
We choose a variational ansatz for the Wigner distribu- 
tion function that is spatially invariant: 
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where 9() is Heaviside's step function. Here the pos- 
itive parameter a represents deformation of Fermi sur- 
face |l4j, the constant k F is the Fermi wave number and 
is related to the number density through n/ = kp/Qii 2 . 
The choice of ([7]) preserves the number density i.e., 
(2ir)- 3 Jd 3 kf{k) =n f . 

The exchange energy can be rewritten as 
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Here we have used the Fourier transform of the dipolar 
potential V dd (q) = (47r/3)eP(3 cos 2 9 q — 1) where 9 q is the 
angle between the momentum q and the dipolar direction 
(i.e., the z-axis) [l5| . We note that the Hartree direct 
term becomes zero for uniform density distribution of 
fermions because the average over the angle 9 q cancels 
out the interaction effect. 

Using the variation ansatz ([7]), the exchange energy 
can be evaluated analytically and is given by 
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where we have defined the "deformation function" : 
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This integral has rather complicated analytical form. It 
is more instructive to plot out the function I(x) which 
we show in Fig. [TJ I(x) is a monotonically decreasing 
function of x, positive for x < 1, passing through zero at 
x = 1 and becomes negative for x > 1. The exchange 
energy ([8]) therefore tends to stretch the Fermi surface 
along the z-axis by taking a — > 0. This however comes 
with the expense of the kinetic energy 
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which favors an isotropic spherical Fermi surface (i.e., 
a = 1). The competition between the two will find an 
optimal value of a in the region a £ (0, 1). The dipo- 
lar interaction therefore, through the Fock exchange en- 
ergy, deforms the Fermi surface of the system. This may 
be regarded as the magnetostriction effect in momentum 
space. 

Inhomogeneous case — Let us now turn to a system of 
N atoms confined in a harmonic trapping potential with 
axial symmetry: 
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We choose a variational Wigner function that has the 
same form as in the homogeneous case, i.e., Eq. ([7]). but 
now the Fermi wave number fcp- is no longer a constant 
and has the following spatial dependence: 
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where ay lo — \jTi/mu and uj = (w 2 ^) 1 / 3 . The vari- 
ational parameters (3 and A represent deformation and 
compression of the dipolar gas in real space, respec- 
tively. Using N = Jd 3 rn(r) one can easily find that 
J:f = (48iV) 1 / 6 A 1 / 2 /a/i . The corresponding density dis- 
tributions in real and momentum space are give by 
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respectively, where R F = (48A r ) 1/6 a,,o/A 1/2 . 

Under this ansatz, each term in the energy functional 
can be evaluated analytically with the total energy given 
by, in units of N 4 / 3 Kui, 
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where c x = 3 1 / 3 /2 8 / 3 ~ 0.2271, c 2 = 2 10 /(3 7 / 2 • 5 • 7tt 2 ) ~ 
0.0634, c dd = d 2 /(a\jiuj), and /3 = (tu r /uj z ) 2/3 mea- 
sures the trap aspect ratio. Here the two terms in the 
square bracket at rhs represent the kinetic and trapping 
energy, respectively, while those in the curly bracket are 
the direct and exchange interaction terms, respectively. 

It is not difficult to see that Eq. fTO]) is not bounded 
from below, a result arising from the fact that the dipolar 
interaction is partially attractive. There however exists, 
under certain conditions, a local minimum in (|10j) . repre- 
senting a metastable state. This situation is reminiscent 
of the case of a trapped attractive BEC [16]. For the 
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FIG. 2: Density distributions in real space (upper plots, in 
units of kp) and in momentum space (lower plots, in units 
of Rp) for an oblate trap with (3q = 0.5 (left plots) and a 
prolate trap with (3q = 2 (right plots). Here r — y 1 x 2 + y' 2 , 
k r = \/k% + kl and c^N 1 ' 6 = 1.5. 



metastable state, the variational energy e(a, /?, A) satis- 
fies the Virial theorem 2Ek in — 2E tr + 3(Ed + E ex ) = 0. 
Hereafter, we refer to the metastable state as the ground 
state. 

We find the ground state by numerically minimizing 
Eq. (fTOf . Fig. [2] shows the ground state density distri- 
butions in both real and momentum space for two dif- 
ferent traps. One can see that while the spatial density 
distributions are essentially determined by the trap ge- 
ometry, the momentum density distributions by contrast 
are quite insensitive to the trapping potential and are in 
both cases elongated along the dipolar direction. Fur- 
ther, the momentum central density at k = decreases 
as (3q increases. 

The stretch in k z is more clearly illustrated in Fig.[3fa), 
where we have plotted the ratio of the root mean square 
momentum in k z direction, yj (k 2 )o, and to that in k x di- 
rection, yj (fc^)o, as a function of the trap aspect ratio (3q 
for several dipolar strengths. It turns out that the dipolar 
interaction leads to nonspherical momentum distribution 
stretched along the dipolar direction irrespective to the 
geometry of trapping potential. This can be attributed 
to the Fock exchange energy that becomes negative for 
a < 1 as discussed in the homogeneous system. This 
result is in stark contrast to the case of dipolar BEC in 
which the Fock exchange energy is absent and the shape 
of the momentum distribution is related to that of the 
spatial distribution through the Fourier transformation. 
Note that, for non-interacting fermions, the resulting mo- 
mentum distribution is isotropic independent of the trap- 
ping potential [13, H] ■ 

Figure 2] shows the real space Thomas- Fermi surface of 
the ground state for different trap geometry. The surface 
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FIG. 3: (Color online) (a) Aspect ratio in momentum space 

v/ {kz)o/ {kx)o and (b) aspect ratio in real space yl {z 2 )o/{x 2 )o 
normalized to that for a non-interacting system as functions 
of /3o for c dd N 1/(i = (solid line), 0.5 (dashed line), 1 (dotted 
line), and 1.5 (dot-dashed line). 
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FIG. 4: (Color online) Thomas-Fermi surface in real space of 
the ground state for noninteracting case (dashed line) and for 
interacting case CddN 1 ^ 6 = 1.5 (solid line). 



of non-interacting fermions has been plotted as dashed 
lines for comparison. While the shape of the cloud re- 
lies on the trap geometry, the dipolar interaction tends 
to stretch the gas along the dipolar direction also in real 
space while compress the gas along the perpendicular ra- 
dial direction. However, once the trapping potential be- 
comes highly elongated (i.e., (3q ^> 1), the dipolar inter- 
action tends to shrink the whole cloud in both the radial 
and the axial directions as shown in the case for (3q = 5 
in Fig. [4j This is because, for such a cigar-shaped trap, 
a number of dipolar fermions align in the axial (dipolar) 
direction and feel strong mutual attractions. 

To better quantify the real space deformation, we show 
in Fig. [3^b) the aspect ratio of the cloud y/ (z 2 ) / (x 2 )o, 
normalized to that of non-interacting Fermi gas, for dif- 
ferent trap geometry. The deviation of the aspect ra- 
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FIG. 5: Critical dipolar interaction strength CddN 1 ^ 6 as a 
function of /3q. 



tio from non-interacting case is most dramatic for near 
spherical traps with 0q w 1. 

As we have already pointed out, only metastable dipo- 
lar gas can exist in a trap. For sufficiently strong dipolar 
interaction strengths, Eq. (JTUJ) no longer supports the lo- 
cal minimum and the system is expected to be unstable 
against collapse. For a given trap aspect ratio /3 , the 
onset of this instability gives rise to a critical value for 
CddN 1 / 6 , yielding the phase diagram as shown in Fig. [5] 
In Ref. [7], Goral et al. showed that, for a sufficiently 
oblate trap, the system is always stable regardless of the 
strength of the dipolar interaction. The critical trap as- 
pect ratio they found is, when translating to our nota- 
tion, about fto « 0.33. Our calculation, however, indi- 
cates that no such critical value exists in f3 — when the 
trap becomes more and more oblate, the critical dipo- 
lar strength increases rapidly but always remains finite. 
This contradiction between our work and theirs origi- 
nates from the Fock exchange term of the dipolar inter- 
action that plays a significant role near critical point and 
is not properly treated in Ref. 0- 

Finally, let us briefly discuss how to detect the dipo- 
lar effects in Fermi gases. In fact, the momentum space 
magnetostriction effect indicates that the dipolar effects 
will manifest in time-of-flight images. Assuming ballis- 
tic expansion after turning off the trapping potential, for 
time t^>l/uj, the aspect ratio of the expanded cloud ap- 
proaches the initial in-trap aspect ratio of the momentum 
distribution as shown in Fig. [3] (a). In other words, re- 
gardless of the initial trap geometry, the expanded cloud 
will eventually become elongated in the dipolar direction. 
In comparison, the expansion of a non-interacting Fermi 
gas is isotropic and results in a spherical cloud in the 
long time limit 



also been used to detect the dipolar effects in chromium 
BEC. In that case, the expansion dynamics strongly de- 
pends on the initial trap geometry 19|, 2fj| - 

In conclusion, we have analyzed the properties of a 
dipolar Fermi gas using a variational method. We have 
emphasized the important roles played by the Fock ex- 
change energy. We found that the dipolar interaction in- 



duces deformation of the Fermi surface and of the phase 
space density distribution. The resulting anisotropic mo- 
mentum distribution of the dipolar gas, which can be 
readily probed using time-of-flight technique, is elongated 
in dipolar direction irrespective to the trap geometry. 

Future work will extend these considerations into col- 
lective excitations and superfluidity of the dipolar gas. 
Since Fermi surface is a key ingredient for low energy 
properties of Fermi gases, a non-spherical Fermi surface 
will lead to new features in collective phenomena in dipo- 
lar Fermi gases. 
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